Enhancer only peaks (well, removing peaks in the TSS)
library(tidyverse)
peak <- read_tsv('/Volumes/data/projects/nei/hufnagel/iPSC_RPE_ATAC_Seq/peak_full/homer_enhancers/all_common_peaks.blackListed.narrowPeak.closestTSS__interesting_homer_motif.bed.gz', col_names = F)
Parsed with column specification:
cols(
X1 = col_character(),
X2 = col_integer(),
X3 = col_integer(),
X4 = col_character(),
X5 = col_integer(),
X6 = col_character(),
X7 = col_integer(),
X8 = col_integer(),
X9 = col_character(),
X10 = col_integer(),
X11 = col_character(),
X12 = col_character(),
X13 = col_character(),
X14 = col_integer(),
X15 = col_integer(),
X16 = col_character(),
X17 = col_double(),
X18 = col_character()
)